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Spin-orbit coupling (SOC), the intrinsic interaction between a particle spin and its motion, is responsible for 
various important phenomena, ranging from atomic fine structure to topological condensed matter physics. 
The recent experimental breakthrough on the realization of SOC for ultra- cold atoms provides a completely 
new platform for exploring spin-orbit coupled superfluid physics. However, the SOC strength in the 
experiment is not tunable. In this report, we propose a scheme for tuning the SOC strength through a fast 
and coherent modulation of the laser intensities. We show that the many-body interaction between atoms, 
together with the tunable SOC, can drive a quantum phase transition (QPT) from spin-balanced to 
spin-polarized ground states in a harmonic trapped Bose-Einstein condensate (BEC), which resembles the 
long-sought Dicke QPT. We characterize the QPT using the periods of collective oscillations of the BEC, 
which show pronounced peaks and damping around the quantum critical point. 

SOC plays a major role in many important condensed matter phenomena and applications, including spin 
and anomalous Hall effects 1 , topological insulators 2 , spintronics 3 , spin quantum computation, etc. In the 
past several decades, there has been tremendous efforts for developing new materials with strong SOC and 
new methods for tuning SOC with high accuracy for spin-based device applications 4 ' 5 . However, the SOC strength 
in typical solid state materials (e.g., ~10 4 m/s in semiconductors) is generally much smaller than the Fermi 
velocity of electrons (~10 6 m/s), and its tunability is also limited and inaccurate. 

On the other hand, the recent experimental breakthrough on the realization of SOC for ultra- cold atoms 6 
provides a completely new platform for exploring SOC physics in both BEC 7-12 and degenerate Fermi gases 13 " 16 . In 
a degenerate Fermi gas, such SOC strength can be at the same order as (or even larger than) the Fermi velocity of 
atoms. Because of the strong SOC, spins are not conserved during their motion and new exotic superfluids may 
emerge. For instance, new ground state phases (e.g., stripes, phase separation, etc.) may be observed in spin-orbit 
coupled BEC 8 " 12 and new topological excitations (e.g., Weyl 14 and Majorana 13 fermions) may appear in spin-orbit 
coupled Fermi gases. The observation and applications of these exciting phenomena require fully tunable SOC for 
cold atoms to characterize the effects of SOC in various phases. Unfortunately, the strength of the SOC in the 
experiment 6 and other theoretical proposals 17 " 20 is not tunable because the SOC strength is determined by the 
directions and wavelengths, not the intensities, of the applied lasers. 

In this report, we propose a scheme for generating tunable SOC for cold atoms through a fast and coherent 
modulation of the Raman laser intensities 21 , which can be easily implemented in experiments. Such tunable SOC 
for cold atoms provides a powerful tool for exploring new exotic Bose and Fermi superfluid phenomena. Here we 
focus on a quantum phase transition (QPT) 22 in a harmonic trapped BEC induced by the many-body interaction 
between atoms and the tunable SOC strength. With the increasing SOC strength, there is a sharp transition for the 
ground state of the BEC from a spin balanced (i.e., equally mixed) phase to a spin fully polarized phase beyond a 
critical SOC strength (i.e., the quantum critical point). By mapping the spin-orbit coupled interacting BEC to the 
well-known quantum Dicke model 23 ' 24 , we obtain analytic expressions for the quantum critical point and 
the corresponding scaling behaviors for the QPT, which agree well with the numerical results obtained from 
the mean-field Gross-Pitaevskii (G-P) equation for the BEC. 

The realization of QPT in the Dicke model using the spin- orbit coupled BEC opens the door for many 
significant applications in quantum optics, quantum information, and nuclear physics 25 " 27 . Previously the 
Dicke model has been studied in several experimental systems 28 " 30 , especially atoms confined in an optical cavity. 
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However the coupling between atoms and optical cavity fields is very 
weak, and the experimental observation of the QPT in the Dicke 
model only occurred recently using the momentum eigenstates for 
a BEC confined in a cavity 31 . Compared with the cavity scheme, the 
spin-orbit coupled BEC utilizes the many-body interaction between 
atoms and has the advantage of essentially no dissipation, fully tun- 
able parameters, very strong coupling, and the use of atom internal 
states, thus provides an excellent platform for exploring Dicke model 
related applications. 

Finally, the QPT is characterized using collective oscillations of the 
BEC, such as the center of mass (COM) motion and the scissors 
mode, where the oscillation periods show pronounced peaks at the 
quantum critical point. Furthermore, the oscillations of the BEC have 
regular periodic patterns in both spin balanced and polarized phases, 
but show strong damping in the transition region. 

Results 

System and hamiltonian. The harmonic trapped BEC in conside- 
ration is similar as that in the recent benchmark experiment 6 . For 
simplicity, we consider a two-dimensional (2D) BEC in the xy plane 
with a strong confine ment (with a trapping frequency co z ) along the z 
direction. Such 2D setup does not affect the essential physics because 
the z direction is not coupled with the SOC. Two hyperfine ground 
states ||)= \F= l y m F = -1) and ||) = \F = l,m F = 0) of 87 Rb atoms 
define the spins of atoms, which are coupled by two Raman lasers 
(with Rabi frequencies ft x and Q 2 ) incident at a n/4 angle from the x 
axis, as illustrated in Figs, la and lb. 

The dynamics of the BEC are governed by the nonlinear G-P 
equation 

ihd<S> /dt = (p 2 /2m + V(r) + Hs + Hj) O, (1) 

under the dressed state basis |f ) = exp(/k 1 -r)| / f), |J)=exp 
(/k 2 -r)| I), where k x and k 2 are the wavevectors of the lasers. <I> = 
(<$)p 0^) T is the wavefunction on the dressed state basis and satisfies 
the normalization condition fdxdy( | ®j | 2 + |®j| 2 ) = 1. The harmonic 
trapping potential V(r) = \ mco 2 (rj 2 x 2 +y 2 ), where co y is the trapping 
frequency in the y direction, and r\ = co x /co y is the ratio of the trap- 
ping frequencies. H s = yp x o z + ti£lo x l2 is the coupling term induced 
by the two Raman lasers with o z and o x as the Pauli matrices. The 
SOC strength y = tried m , k L = \Ki -K 2 |/2 = n / 1, and y is the 
wavelength of the Raman lasers. The Raman coupling constant 
Q = QiQl/2A with A as the detuning from the excited state. The 
mean field nonlinear interaction term Hj = diag (g^ | <X> T | 2 + 1 ® j | 2 , 
gjj|<I>j| 2 + gjjl^jl 2 ), where the inter-and intraspin interaction con- 
stants gjj = gn = 47ifj 2 N(c 0 + c 2 )lma z and g^ = 47zft 2 Nc 0 /ma zy c 0 and 
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Figure 1 | An illustration of the experimental scheme for realizing SOC 
for cold atoms 6 , (a) Laser setup, (b) Atom laser coupling, (c) A typical 
density distribution of one spin component of the BEC in the ground state, 
(d) Corresponding momentum distribution. The vertical and horizontal 
dotted lines are k x = 0 and k y = 0 respectively. 



c 2 describe the corresponding s-wave scattering lengths 32 , N is the 
atom number, and a z = y / 2nh/mco z . 

Because the SOC strength y is determined by the laser wavevector 
k L , the SOC energy can be comparable to or even larger than other 
energy scales (e.g., the Raman coupling Q) in the BEC. In a Fermi gas, 
y can be larger than the Fermi velocity of atoms. Unfortunately, due 
to the same reason, y cannot be easily adjusted in experiments, which 
significantly restricts the applications of the SOC in cold atoms. Note 
that although theoretically it may be possible to tune the SOC 
strength by varying the angle between two Raman beams, experi- 
mentally it is impractical because of many limitations of the experi- 
mental setup. 

Tunable SOC for cold atoms. We propose a scheme for tuning the 
SOC strength y through a_fast and coherent modulation of the 
Raman coupling ft = ft 0 + ft cos(a)£) that can be easily realized in 
experiments by varying the Raman laser intensities. For ft >ft 0 , ft 
changes sign at certain time, which can be achieved by applying a n 
phase shift on one Raman laser. Here the modulation frequency co is 
chosen to be much larger than other energy scales in Eq. (1). In 
this case, the Hamiltonian in Eq. (1) can be transformed to a 
time-independent one using a unitary transformation \j/ = exp 
[/Qsin(o)f)(T x /(2ct})]0. After a straightforward calculation with the 
elimination of the fast time-varying part in the Hamiltonian 33 ' 34 , the 
nonlinear G-P equation (1) becomes 

ihdi///dt= [p 2 /2tn + V(r)+H s + H I ]xl/, (2) 

where the Raman coupling becomes H s = y e ffp x (J z + fiQQG x /2 with 
the effective SOC strength 

y eff = yj 0 (h/w). (3) 

Here / 0 is the zero order Bessel function. Clearly, y e jj can be tuned 
from the maximum y without the modulation to zero with a strong 

modulation. The mean field interaction term Hj = a ^ | | 2 + | xjj^ | 2 ^ 

\jf + prijj, a = j, p = (gn — f )/2 and T is a 2 X 2 matrix whose 

i \2 [3 

elements are given by rn = |i/^| 
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Quantum phase transition. The tunable SOC, in combination with 
the many-body interaction between atoms, can drive a quantum 
phase transition between different quantum ground states in a 
harmonic trapped BEC. Here the ground state of the BEC is 
obtained numerically through an imaginary time evolution of the 
G-P equation (2). A typical density profile of the ground state is 
shown in Fig. lc, which has a Thomas-Fermi shape, similar as that 
in a regular BEC. However, the momentum distribution of the BEC 
has a peak around the single particle potential mini mum located at 
Ky) = (-K min , 0) (see Fig. Id), where K min = ^y 2 eff -Q 2 0 /4y 2 eff 

and the degeneracy between ±iC min is spontaneously broken. 

To characterize the ground state of the spin- orbit coupled BEC, we 
calculate the spin polarization \{o z )\ = \f dr(\i//^\ 2 — |iAJ 2 )|, an d 
(a x ) =2Refdri//t\l/^. Here we choose the absolute value of (a z ) 
because the two degenerate ground states at ±K min have opposite 
(a z ) due to the spinmomentum locking termp^ and they are spon- 
taneously chosen in experiments. In Fig. 2a, we plot \(<J Z )\ and {o x ) 
with respect to y e ff. For a small y e p the spin up and down atoms have 
an equal population, thus (a z ) = 0, (a x ) = — I, i.e., the spin balanced 
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Figure 2 | Quantum phase transition with tunable y eff . Q 0 = 16, y = 9.37 
is the bare SOC strength without modulation, (a) Plot of the spin 
polarization I (er z ) I and (cr x ) in the ground state. The blue and red lines are 
from the prediction of the Dicke Hamiltonian. The circles and squares are 
from the numerical simulation of the G-P equation (2). (b) Plot of the 
COM motion period T. The shift of the harmonic trap D = 1, T 0 = 2%lco y . 

phase. Beyond a critical point y^-, the spin population imbalance 
increases dramatically and reaches the spin polarized phase (a z ) = 
1 = 0) within a small range of y e ff. The spin balanced and spin 
polarized phases at small and large y e jj can be understood from the 
single particle Hamiltonian H s , where the Raman coupling Q 0 (r x /2 
and SOC y e jjp x (J z dominate at the small and large y e ff respectively. 
More nu merical results show that the critical transition point occurs 

aty< =y / £V2. 

The QPT from spin-balanced to spin polarized phases can be 
understood by mapping the spin-orbit coupled BEC to the quantum 
Dicke model. For an interacting BEC in a harmonic trap with a large 
atom number N (so that the mean field theory works), all atoms are 
forced to occupy the same many-body ground state (i.e., the state in 
Fig. lc). Therefore the energy variation for the change of the spin 
(e.g., spin flip) of one atom need be determined by the coupling 
between the atom spin and the many-body ground state mode. 
This is very different from an non-interacting BEC where atoms 
do not affect each other, but the same as that for many atoms inter- 
acting with a single photon mode in an optical cavity 23 . Treating the 
interacting many-body ground state as a single mode composed of 
different harmonic trap modes, we can map the Hamiltonian for the 
spin-orbit coupled BEC to 

, y pff \/mhoj x , , 

H = Pko x N^ a + m 0 S x + m (a 1 " - a) (S + - S _ ) , (4) 

V2 

which is similar to the Dicke model for two -level atoms coupled with 
a cavity field 23 . S x , yyZ are the large spins for all atoms, S+ = S y + iS z , S- 
= S y — iS z , a f a is a harmonic trap mode, a = y / mco x /2h 
(x + ip x /m(D x ). In the mapping to the Dicke Hamiltonian (4), we 
neglect two minor effects of the interaction terms in Eq. (2): 1) a 
constant mean field background energy term that does not affect the 
dynamics; 2) a small term kS 2 z /N 2 with k ~ g^ — g^. The second 
term can shift the Dicke phase transition point. However, for realistic 
experimental parameters for 87 Rb, k is very small and the effects of 
kS 2 z /N 2 can be neglected, as we confirm in the numerical simulation 
of the G-P equation. Nevertheless, these minor effects cannot be 
captured by the non-interacting BEC. The QPT from spin-balanced 
to spin-polarized phases in the spin-orbit coupled BEC is similar to 



that from normal to superradient phases in the Dicke model. The 
critical point for the QPT can be derived from the standard mean- 
field approximation 24 , yielding the relation y c e ^ = ^/Q 0 /2, which is 
exactly the same as that from numerically simulating the G-P equa- 
tion (2). Clearly, the QPT can also be driven by varying Q 0 for a fixed 
y e ff. Just beyond the critical point y c e p the Dicke model predicts that 
the scaling of the order parameters is \{g z )\=2\S z \/N = 

\J l -(y C e ff /yefft> (°x)=2S x /N=-(y c eff /y ejf ) 2 for y eff >y c eff , 
and \{<J Z )\ = 0, (<j x ) = —1 for y e ^<y c e ^. Such scaling behaviors are 
confirmed in our numerical simulation of the G-P equation (see 
Fig. 2a). The perfect match between numerical results from the G- 
P equation and the predictions of the Dicke Hamiltonian shows the 
validity of the mapping to the Dicke model. 

We emphasize that the many-body interaction between atoms in 
Eq. (2) plays a critical role in the QPT by forcing all atoms in a single 
spatial mode. For a non-interacting BEC, the atoms can occupy both 
— ^min in the momentum space with an artificial ratio because these 
two states are energetically degenerate and there is no correlation 
between atoms. The resulting spatial distribution of the BEC is thus 
artificial and the above single spatial mode approximation in Eq. (4) 
does not apply. Our numerical simulation of the G-P equation with- 
out interactions also shows {o z ) = 0 or other random values in 
certain region of y e ^>y c e ^, which disagrees with the prediction of 
the Dicke model. This disagreement confirms that atoms in a non- 
interacting BEC do not response to the change of y e ff collectively, 
although non-interacting and interacting BECs share the same trans- 
ition for the energy spectrum at y c e ^, which changes from one single 
minimum at K x = 0 to two minima at ±K min . While for interacting 
BECs with large atom numbers N = 4 X 10 4 and 10 6 , we obtain 
exactly the same results as that in Fig. 2a, which further confirm 
the validity of our mapping to the Dicke model in the large N limit. 

Collective dynamics in BEC: the signature of QPT. It is well-known 
that various physical quantities may change dramatically around the 
quantum critical point (i.e., critical phenomena), which provides 
additional experimental signatures of the QPT. We focus on two 
types of collective dynamics of the ground state of the BEC: the 
COM motion and the scissors mode induced by a sudden shift or 
rotation of the harmonic trapping potential, respectively. In a regular 
BEC without SOC, the COM motion is a standard method to 
calibrate the harmonic trapping frequency because the oscillation 
period depends only on the trapping frequency 35 and is not 
affected by other parameters such as nonlinearity, shift direction 
and distance, etc. 

We numerically integrate the G-P equation (2) and calculate the 
COM <r(f)> = jdxdy(\\jj^(t)\ 2 + |iA|M| 2 )rW. The COM motion 
strongly depends on the direction of the shift D of the harmonic 
trap. When D is along they direction, the period of the COM motion 
along they direction is T 0 = 2nl(D y and not affected by y e ff, while the 
COM motion in the x direction disappears (i.e., (x) = 0). Here the 
COM period Tis obtained through the Fourier analysis of (r (t)). The 
physics is very different when D is along the x direction, where (y (£)) 
= 0 as expected, but {x (t)) depends strongly on y e p as shown in 
Fig. 2b. In Fig. 3, we also plot T as a function of y^and Q 0 . Without 
SOC (y e ff = 0), T = T 0 , the period for a regular BEC, as expected. T 
increases with y e jj in the spin-balanced phase, but decreases when 
spin starts to be polarized, leading to a sharp peak at the quantum 
critical point y^. The oscillation of (x(t)) in the spin balanced phase 
is completely clissipationless, while a strong damping occurs in a 
small range of y e ff beyond y^ (see the inset in Fig. 3). Far beyond 
y^, the oscillation becomes regular again with the period T = T 0 
because the ground state has only one component in this region. The 
peak and the damping of the oscillation around y^ provide clear 
experimental signatures for the QPT. Moreover, T also depends on 
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Figure 3 | Plot of the COM motion period T versus / e //and fl 0 . The insets 
show the corresponding (x ( t)). The green line is the theoretical prediction 
y c e ff = a/Q 0 /2 from the Dicke model. The circles are the numerical results 
from the G-P equation (2). T 0 = 2n/w yy D = 1. 

the magnitude D of the shift near the critical point y c e ^: the larger D, 
the smaller T. 

Another collective dynamics, the scissors mode 36 , shows a similar 
feature as the COM motion. The scissors mode can be excited by a 
sudden rotation of the asymmetric trapping potential (i.e., y\¥* 1) by 
an angle 0, which induces an oscillation of the quantity (xy) = 
j dxdy{ | (t) | 2 + (t)\ 2 xy. Without SOC, the period of the scissors 

mode is T\ = 2n j ^Jcol + coj , as observed in experiments 37 . In 

Fig. 4, We plot the oscillation period T with respect to y e jj for three 
different Q 0 - We have confirmed that the same QPT occurs for \(a z )\ 
and (a x ) of the ground state in this as ymmetric potential with the 
quantum critical point y c e ^ = y/Qo/2, as predicted by the Dicke 
model. Similar as the COM motion, we observe the peak and damp- 
ing of the oscillation around y c e ^. Far beyond y^, the oscillation 
period is 7\. Similar as the dependence of the COM motion on the 
shift distance D, the angle 0 also influences the period of the scissors 
mode near y^: the smaller 0, the larger T. 

Discussion 

In summary, we show that the SOC strength in the recent break- 
through experiment for realizing SOC for cold atoms can be tuned 
through a fast and coherent modulation of the applied laser intens- 
ities. Such tunable SOC provides a powerful tool for exploring 




Figure 4 | Plot of the scissors mode oscillation period Twith respect to / e #- 
for three different fl 0 . T\ = 2n j w y ^/r} 2 + l is the oscillation period 
without SOC, rj = y/5 y 6 = 4°. The inset is a typical density distribution of 
the ground state of the BEC. The dotted lines show the corresponding 
quantum critical points predicted by the Dick model. 



spin-orbit coupled superfluid physics in future experiments. By vary- 
ing the SOC strength, the many-body interaction between atoms can 
drive a QPT from spin balanced to spin polarized ground states in a 
harmonic trapped BEC, which realizes the long-sought QPT from 
normal to superradient phases in the quantum Dicke model and may 
have important applications in quantum information and quantum 
optics. 

Methods 

We choose the physical parameters to be similar as those in the experiment 6 : (co y ,co z ) 
= In X (40, 400) Hz, r\ = 1, X = 804.1 nm, c 0 = 100.86a B , c 2 = ~0A6a B 38 with the 
Bohr radius a B ,N = 1 X 10 4 , co = 2n X 4.5 kHz. For the numerical simulation, we 
need a dimensionless G-P equati on that is o btained by choosing the units of the 
energy, length and time as hco y , Jft/ (mcOy) = 1.7 fim, and \lco y — 4 ms, respectively. 
The dimensionless parameters in the G-P equation become y = Jft/ (mco y )k L = 9.37, 

a = 2Ny/2nmco z /h(c 0 + c 2 ) = 495 and />= — Ny/2nmo) z /hc 2 = 1.14. 

Note added. After our manuscript was initially posted at arXiv (arXiv: 1111 .4778), our 
proposed peaks of the dipole oscillation periods (Figs. 2b and 3) were observed 
experimentally 39 , and more detailed theoretical studies were also performed 40 . 

1. Xiao, D., Chang, M.-C. &Niu, Q. Berry phase effects on electronic properties. Rev. 
Mod. Phys. 82, 1959 (2010). 

2. Hasan, M. & Kane, C. L. Topological insulators. Rev. Mod. Phys. 82, 3045 (2010). 

3. Zutic, I., Fabian, J. & Das Sarma, S. Spintronics: fundamentals and applications. 
Rev. Mod. Phys. 76, 323 (2004). 

4. Nadj-Perge, S., Frolov,S. M., bibinfoauthorBakkers, E. P. A. M. & Kouwenhoven, 
L. P. Spin-orbit qubit in a semiconductor nanowire. Nature 468, 1084 (2010). 

5. Jespersen, T. S. et al. Gate-dependent spinorbit coupling in multielectron carbon 
nanotubes. Nature Phys. 7, 348 (2011). 

6. Lin, Y.-J., Jimenez- Garcia, K. & Spielman,I. B. Spin-orbit-coupled Bose-Einstein 
condensates. Nature 471, 83 (2011). 

7. Wu, C., Mondragon-Shem, I. & Zhou,X.-F. Unconventional Bose-Einstein 
condensations from spin-orbit coupling. Chin. Phys. Lett. 28, 097102 (2011). 

8. Wang, C., Jian, C.-M. & Zhai, H. Spin-orbit coupled spinor Bose-Einstein 
condensates. Phys. Rev. Lett. 105, 160403 (2010). 

9. Ho, T.-L. & Zhang, S. Bose-Einstein condensates with spin-orbit interaction. Phys. 
Rev. Lett. 107, 150403 (2011). 

10. Zhang, Y., Mao, L. & Zhang, C. Mean-field dynamics of spin-orbit coupled Bose- 
Einstein condensates. Phys. Rev. Lett. 108, 035302 (2012). 

11. Hu, H., Ramachandhran, R, Pu, H. & Liu, X.-J. Spin-orbit coupled weakly 
interacting Bose-Einstein condensates in harmonic traps. Phys. Rev. Lett. 108, 
010402 (2012). 

12. Sinha, S., Nath, R. & Santos, L. Trapped two-dimensional condensates with 
synthetic spin-orbit coupling. Phys. Rev. Lett. 107, 270401 (2011). 

13. Zhang, C., Tewari, S. Lutchyn, R. M. & Das Sarma, S. p x + ip y Superfluid from s- 
wave interactions of Fermionic cold atoms. Phys. Rev. Lett. 101, 160401 (2008). 

14. Gong, M., Tewari, S. & Zhang, C. BCS-BEC crossover and topological phase 
transition in 3D spin-orbit coupled degenerate Fermi gases. Phys. Rev. Lett. 107, 
195303 (2011). 

15. Hu, H., Jiang, L., Pu, H. & Liu, X.-J. Probing anisotropic superfluidity in atomic 
Fermi gases with Rashba spin-orbit coupling. Phys. Rev. Lett. 107, 195304 (201 1). 

16. Yu, Z.-Q. & Zhai, H. Spin-orbit coupled Fermi gases across a Feshbach resonance. 
Phys. Rev. Lett. 107, 195305 (2011). 

17. Ruseckas, J., Juzeliunas, G., Ohberg, P. & Fleischhauer, M. Non-Abelian gauge 
potentials for ultracold atoms with degenerate dark states. Phys. Rev. Lett. 95, 
010404 (2005). 

18. Zhang, C. Spin-orbit coupling and perpendicular Zeeman field for fermionic cold 
atoms: observation of the intrinsic anomalous Hall effect. Phys. Rev. A 82, 
021607(R) (2010). 

19. Dalibard, J., Gerbier, F., Juzeliunas, G. & Ohberg, P. Artificial gauge potentials for 
neutral atom. Rev. Mod. Phys. 83, 1523 (2011). 

20. Campbell, D. L., Juzeliunas, G. & Spielman, I. B. Realistic Rashba and Dresselhaus 
spin-orbit coupling for neutral atoms. Phys. Rev. A 84, 025602 (2011). 

21. Grifoni, M. & Hanggi, P. Driven quantum tunneling. Phys. Rep. 304, 229 (1997). 

22. Sachdev, S. Quantum Phase transtitions. (Cambridge University Press, 
Cambridge, England, 1999). 

23. Dicke, R. H. Coherence in Spontaneous Radiation Processes. Phys. Rev. 93, 99 
(1954). 

24. Emary, C. & Brandes, T. Quantum chaos triggered by precursors of a quantum 
phase transition: the Dicke model. Phys. Rev. Lett. 90, 044101 (2003). 

25. Gross, M. & Haroche, S. Superradiance: an essay on the theory of collective 
spontaneous emission. Phys. Rep. 93, 301 (1982). 

26. Klein, A. & Marshalek, E. R. Boson realizations of Lie algebras with applications to 
nuclear physics. Rev. Mod. Phys. 63, 375 (1991). 

27. Hammerer, K., Sorensen, A. S. & Polzik, E. S. Quantum interface between light 
and atomic ensembles. Rev. Mod. Phys. 82, 1041 (2010). 



SCIENTIFIC REPORTS | 3 : 1937 | DOI: 1 0.1 038/srep01 937 



4 



28. Schoelkopf, R. J. & Girvin, S. M. Wiring up quantum systems. Nature 451, 664 
(2008). 

29. Keller, M., Lange, B., Hayasaka, H., Lange, W. & Walther, H. Continuous 
generation of single photons with controlled waveform in an ion-trap cavity 
system. Nature 431, 1075 (2004). 

30. Garraway, B. M. The Dicke model in quantum optics: Dicke model revisited. Phil. 
Trans. R. Soc. A 369, 1137 (2011). 

31. Baumann, K., Guerlin, C., Brennecke, F. & Esslinger, T. Dicke quantum phase 
transition with a superfluid gas in an optical cavity. Nature 464, 1301 (2010). 

32. Ho, T.-L. Spinor Bose condensates in optical traps. Phys. Rev. Lett. 81, 742 (1998). 

33. Eckardt, A., Weiss, C. & Holthaus, M. Superfluid- insulator transition in a 
periodically driven optical lattice. Phys. Rev. Lett. 95, 260404 (2005). 

34. Lignier, H. et al. Dynamical Control of Matter- Wave Tunneling in Periodic 
Potentials. Phys. Rev. Lett. 99, 220403 (2007). 

35. Stringari, S. Collective Excitations of a Trapped Bose-Condensed Gas. Phys. Rev. 
Lett. 77, 2360 (1996). 

36. Guery-Odelin, D. & Stringari, S. Scissors Mode and Superfluidity of a Trapped 
Bose-Einstein Condensed Gas. Phys. Rev. Lett. 83, 4452 (1999). 

37. OMarago, O. M. et al. Observation of the Scissors Mode and Evidence for 
Superfluidity of a Trapped Bose-Einstein Condensed Gas. Phys. Rev. Lett. 84, 2056 
(2000). 

38. Widera, A. et al. Precision measurement of spin-dependent interaction strengths 
for spin-1 and spin-2 87 Rb atoms. New J. Phys. 8, 152 (2006). 

39. Zhang, J.-Y. et al. Collective Dipole Oscillations of a Spin- Orbit Coupled Bose- 
Einstein Condensate. Phys. Rev. Lett. 109, 115301 (2012). 



40. Li, Y., Martone, G. I. & Stringari, S. Sum rules, dipole oscillation and spin 
polarizability of a spin-orbit coupled quantum gas. EuroPhys. Lett. 99, 56008 
(2012). 



Acknowledgements 

We thank helpful discussion with Peter Engels, Li Mao and Chunlei Qu. This work is 
supported by DARPA-YFA (N66001- 10- 1-4025), ARO (W911NF-12-1-0334), and 
NSF-PHY (1104546). Gang Chen is also supported by the 973 program under Grant No. 
2012CB921603 and the NNSFC under Grant No. 11074154 and 61275211. 

Author contributions 

Y.Z., G.C. and C.Z. conceived the idea, Y.Z., G.C., C.Z. performed the calculation, Y.Z., G.C., 
C.Z. wrote the manuscript, C.Z. supervised the whole research project. 

Additional information 

Competing financial interests: The authors declare no competing financial interests. 
License: This work is licensed under a Creative Commons 

Attribution-NonCommercial-ShareAlike 3.0 Unported License. To view a copy of this 
license, visit http://creativecommons.Org/licenses/by-nc-sa/3.0/ 

How to cite this article: Zhang, Y., Chen, G. & Zhang, C. Tunable Spin-orbit Coupling and 
Quantum Phase Transition in a Trapped Bose-Einstein Condensate. Sci. Rep. 3, 1937; 
DOI:10.1038/srep01937 (2013). 



SCIENTIFIC REPORTS | 3 : 1937 | DOI: 1 0.1 038/srep01 937 



5 



